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We show that, despite appearances, a theoretical approach to neutrino oscillation in which 
CIh' the neutrino and its interaction partners are entangled yields the standard result for the 

Oh' neutrino oscillation wavelength. We also shed some light on the question of why plane-wave 

, approaches to the neutrino oscillation problem can yield the correct oscillation wavelength 

even though they do not explicitly account for the localization of the neutrino source and 

CN ■ 

' the detector. 
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\D ■ I. INTRODUCTION 

O 
O 

Is the standard expression jjj] for the wavelength of neutrino oscillation correct? Does this 
' expression depend correctly on the underlying neutrino mass splittings, so that the splittings that 
^ , have been inferred from data using this expression are right? 

The quantum-mechanical physics of neutrino oscillation has proved to be quite subtle, and has 
been analyzed in a variety of ways over the years. There have been treatments based on plane 



waves 



^, on neutrino wave packets [5|-l9|], and on quantum field theory 1CH14||. (A thorough 
summary has been given recently by Akhmedov and Smirnov [3].) Most of these treatments have 
yielded the standard expression for the probability of oscillation, but the correct way to derive 
this expression is still occasionally disputed. The analysis presented in the present work, based 
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on plane waves and on an assumed entanglement between the oscillating neutrino and a recoil 
particle, finds an oscillation wavelength that appears to be very different from the usual one, but 
turns out to be physically equivalent to it. 

In this paper, we recap the recent treatment of oscillation using entangled plane waves. We 
express a number of concerns about such an analysis, and about plane wave treatments in general. 
However, we also comment on why plane wave treatments, while not strictly correct and consis- 
tent, nevertheless often yield the correct oscillation wavelength. We demonstrate that, in spite of 
concerns, and in spite of the fact that the oscillation wavelength found by considering entangled 
plane waves appears to be markedly different from the standard one, this wavelength, properly 
understood, is physically completely equivalent to the usual one, and yields the same neutrino 
mass splittings when applied to data. This result is the main point of this paper. 



II. PLANE- WAVE DESCRIPTION WITH ENTANGLEMENT 



In this section, we present a derivation of the oscillation wavelength, focusing on neutrinos 
produced in two-body decays of the form P — )• + R, in which a parent P decays into a neutrino 
V and a recoil R. In an effort to avoid apparent paradoxes that can arise in standard plane- 
wave treatments of neutrino oscillations, we assume that the neutrino and its recoil form a two- 
body entangled state. The role of entanglement in neutrino and neutral kaon oscillations has 



been considered by others 16H23t|. In the analysis presented here, the oscillation wavelength that 
emerges corresponds to the separation of the neutrino and its associated recoiling particle at a 
common time in the parent rest frame. This wavelength is not the usual wavelength that would 
be observed in an experiment (because the recoiling particle is not usually observed), but we 
show that the two wavelengths are physically equivalent and lead to the same neutrino survival or 
appearance probability.^ 

A neutrino of flavor a is a linear combination of mass eigenstates Vi with masses m,, 

W) = Y.KM). (1) 



where the Uai are elements of the Maki-Nakagawa-Sakata-Pontecorvo mixing matrix In the 
approach that assumes entanglement, energy-momentum conservation requires that the energy 



^ This section incorporates arguments presented previously in unpublished form by one of us (RGHR) [ij]. A 
revised version of that work [25| takes into account insights explained in the present paper. 
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and momentum of the recoil in the decay P ^ z/ + R depend on which Vi was actually emitted. 
In the rest frame of the parent, pn = —pu, where pji and pi, are the momenta of the recoil and 
neutrino, respectively. In terms of the coordinates in that frame, the two-particle wave function 
may be written as 



Here, \R) describes the internal degrees of freedom of the recoil, and simply identifies one of 
the three neutrino mass eigenstates. Ej^^i and Eji^i are the energies of the neutrino and recoil, 
respectively, and t is the time. The P-rest-frame distance between the neutrino and the recoil is 
D = x—X, where the coordinates x and X are the positions of the neutrino and recoil, respectively. 
The energy Ea = E^^i + Ep^^i = mp in the third line of Eq. [2]is just the mass of the parent P. The 
index i on p^j^i, Ey^i, and Eji^i indicates that these quantities all depend on i, while Ea, as well as 
t, X, and X, do not. (To describe a neutrino oscillation experiment with a spatial resolution much 
better than the oscillation length, we have to fix (t, x) and (t, X) at the P-rest-frame coordinates 
of the spacetime points where the neutrino and recoil, respectively, are detected. These points are 
defined by the experiment and therefore do not depend on i.) Note that the third line of Eq. [2] 
shows that, in the rest frame of the parent, the two-particle system consisting of the neutrino 
and the recoil can be described by a one-particle wave function depending only on the relative 
coordinate D. This is analogous to the treatment of the hydrogen atom in most quantum mechanics 
textbooks, where the two-body wave function of the electron and the nucleus is factorized into a 
one-particle wave function describing the relative motion of the two particles, and a one-particle 
wave function describing the center of mass motion (which is trivial in the center of mass frame) . 

Suppose that in the laboratory frame, the parent P is moving to the right along the x-axis 
with speed /3. The lab-frame location x' of the neutrino when it is at the spacetime point (t, x) is 
x' = 7(x + (3t). The corresponding lab-frame location X' of the recoil at the same time t in the 
P rest frame is X' = j{X + fit). Thus, the lab- frame distance between the neutrino and recoil at 
time t in the P rest frame is x' — X' = j{x — X) = D' . Consequently the state \Ri^; x, X) of Eq. [2] 
may be rewritten as 



Rv;x,X) ~ Ul^\Ri{pp^i, X))\vi{pu^ix)) 





(2) 




(3) 
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where p^^i is the parent rest-frame momentum of the neutrino when the latter is fj. The leading 
energy-dependent phase factor in Eq. [3] is unobservable. As a result, the interference effects of 
neutrino oscillations arise solely from the different momenta in the components in the final state. 
In the oscillation probability at lab-frame separation D', P{i'a — ^ i^/?) = {Ri^plRi^a) £)/ , the i — j 
interference term depends on D' through the phase factor exp [i7^^(j)jy_j — pj^j)!)']. Thus, the 
wavelength A'^, of oscillation in the recoil-neutrino separation D' is determined by 

l~\p^,^-P^,j)>^'D,ij = 2vr. (4) 
Now, to leading order in Amfj = mf — m|, 

2 I ^2 



2 1 mp + mr, 

Pu,i-Pu,j = -Arriij- 2 2"' ^ 

■' 2mp nip — mj^ 

where niR is the mass of the recoil. Thus, apart from an irrelevant sign, 

Ai..- - 'rr^^r^l - (6) 

' Am^j nip + m|j 

To compare this result with the standard expression for the wavelength of neutrino oscillation, 
it is useful to rewrite it in terms of the "neutrino beam energy" Eq, defined as the energy that 
massless neutrinos would have in the laboratory. Since Eq = 7(1 + /3)(mp — m^)/2mp, and the 
speed of the recoil in the P rest frame obeys 1 + /3/j = 2m'^p/{w?p + nip), we have 

_ At^E'q I + Pr 

Equations [6] and [7| define a wavelength in laboratory coordinates for the separation D' between 
the neutrino and the recoil. In the standard expression [l| for neutrino oscillation, the survival 
probability oscillates as a function of the lab-frame distance L' between the neutrino source and 
the detector with a wavelength X'^ - given by 

'■J 

We will show in the next section that these two expressions, Eqs. [7]and[8l in fact give equivalent 
results for the laboratory wavelength of neutrino oscillations measured in the standard fashion by 
detecting neutrinos at a known distance from a source region. 

A comment is in order about the role that the entanglement between the neutrino and the 
recoil plays for the oscillation phenomenology. In particular, one may wonder if the oscillation 
pattern is modified if the recoil undergoes an interaction that breaks the entanglement long before 
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the neutrino is detected. In fact, in a typical experiment, the recoil will interact with the matter 
that makes up the neutrino source very soon after it has been produced. If it is unstable, it might 
also decay very rapidly. However, it is easy to see that such interaction or decay cannot change the 
neutrino oscillation phenomenology. First, the order in which the recoil interaction or decay and 
the neutrino detection occur depends on the Lorentz frame in which we are working. The neutrino 
flavor transition probabilities, on the other hand, are Lorentz invariant, so they cannot depend on 
the time ordering of these processes. In a more formal way, the same conclusion can be reached 
by using the quantum amplitude approach to particle mixing probabilities developed in |2(|. This 
approach, which takes the entanglement between the neutrino and the recoil into account and 
confirms the validity of Eq. [2l yields the joint probability for the recoil to interact (or decay) at 
one spacetime point, and the neutrino to produce a charged lepton of a given flavor a at another 
spacetime point. To determine the consequences of this approach for neutrino oscillations, one 
must integrate this joint probability over all possible interaction points of the recoil. But it turns 
out that the joint probability depends on the recoil interaction point (T, X) only through a phase 
factor Q~'^^R,i^~^'Pv,iX ^ analogous to the recoil phase factor we have encountered in the second line 
of Eq. [2j It is easy to show that this factor is independent of which neutrino mass eigenstate has 
been emitted together with the recoil 0, [2^. Thus, it is an overall phase that does not influence 
the quantum interference between neutrino mass eigenstates that leads to neutrino oscillations. 
Hence, the oscillation pattern is independent of where or how soon the recoil interacts or decays, 
destroying its entanglement with the neutrino. 



III. PHYSICAL EQUIVALENCE OF WAVELENGTHS 

Since the wavelengths and A'j;^ in Eqs. [7] and [8] describe oscillations in two different 

variables, D' and L', let us see how these variables are related. We take the decay P — > + R to 
occur at the spacetime point (t, x) = {t',x') = (0,0). Suppose that the neutrino is then detected 
at a subsequent time t in the P rest-frame. Then, in that frame, if the neutrino's speed is Pi,, it 
will have traveled a distance L = j3yt to its point of detection. The P rest-frame separation D 
between the neutrino and the recoil when the neutrino is detected will be -D = Pyt + Put. Thus, 
since the neutrino is ultrarelativistic, with /3,y ~ 1, 

^ = 1 + /3r- (9) 
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(11) 



Now, in the laboratory frame, the distance L' between the neutrino source (i.e., the location of P 
when the decay P — > + R occurred) and the point of detection is L' = ^{L + fit) ~ 7(1 + /3)L, 
where we have once again used ~ 1. Furthermore, as we have already seen, at the neutrino 
detection time t in the P rest-frame, the laboratory-frame distance D' between the neutrino and 
the recoil is related to its P rest-frame counterpart D hy D' = x' — X' = j{x — X) = 'jD. Thus, 
the variables D' and L' to which the wavelengths A'^, and X'^ refer are related by 

L' 7(1 + (3)L 1 + /3 ^ ^ 

By comparison, from Eqs. [7] and El 

A^ ^ 1 + /3r 

That is, the ratio between the "new" wavelength A^ - of oscillation in D', and the standard 
wavelength A^ of oscillation in L' , is exactly the same as that between D' and L' . 

Thus, although the new A^ of Eq. [7]has emerged from an approach that entails entanglement, 
while the standard A^ of Eq. [8] has come from analyses that generally do not, these two wave- 
lengths differ only because they refer to two different, alternative distance variables. Thus, they 
are physically equivalent. Properly used to fit given oscillation data, they would yield precisely 
the same neutrino squared-mass splitting Am?^-. 

In practice, of course, an analysis of experimental data using Eq. [7| instead of the usual Eq. [8] 
would require knowledge of the coordinates of both the neutrino and the recoil at the same time 
t in the parent rest frame. The required coordinates of the recoil would be extremely difficult to 
obtain. However, it is not necessary to actually detect the recoil if the spacetime trajectory and 
energy of the parent particle are known — as they would be, for instance, in an electron-capture 
beta-beam neutrino experiment. In that case, energy-momentum conservation could be used to 
infer the trajectory of the recoil from the kinematics of the parent and the measured coordinates 
of the neutrino. 



IV. WHY IS THE SAME Am^ OBTAINED? 



So long as the standard approaches and the approach that invokes entanglement are all valid 
ways of deriving the neutrino oscillation wavelength, the wavelengths derived by these approaches 
obviously must be correct and physically equivalent. In the previous section, we have demonstrated 
the physical equivalence for a particular approach that invokes entanglement and that obtains the 
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wavelength for oscillation in the separation D' between the neutrino and its recoiling partner at a 
given time t in their common parent's rest frame. But this must be a more general result. Consider 
a "gedanken" neutrino oscillation experiment, with some relevant properly defined distance D" 
between the neutrino and its recoiling partners, and with D" proportional to the usual source-to- 
detector distance L' . Then, the ratio between the wavelengths of oscillation in D" and L' will be 
D" /L' , so that these two wavelengths will be physically equivalent. 

Another way of understanding intuitively why the standard expression for the oscillation length 
is obtained even when the neutrino and the recoiling particle are considered as an entangled state 
is the following: Let us choose a frame in which the component of the recoil that is entangled with 
neutrino mass eigenstate of mass m\ is at rest; i.e. i = 0. We denote kinematic quantities in 
this frame by a double prime. Prom simple kinematic arguments, it follows that 

Pkj = 2' 2 + • • ■ (12) 
'■^ rrip — m|j 

E'kj = mR + ..., (13) 



where '. . . ' denotes terms that are at least 4-th order in the neutrino masses, divided by combi- 
nations of rriR and mp. The complex phase of the jth component of the entangled state has the 
form 

- E'k/' + p'ijX" - E'l/' + . (14) 

We see that, at leading order in Am|^, the first term is a constant that does not contribute to the 
phase differences relevant to neutrino oscillations and can therefore be omitted. We will now show 
that the second term can be neglected as well. Since we are using a plane wave approach here, 
it is not immediately obvious what X" is; after all, plane waves are delocalized over space. We 
will argue below that it is reasonable to impose the relations x" = v'^t" and X" = v'^t" by hand, 
where it is justified to take v'l and v'^ to be the averages of the group velocities associated with the 
individual components of the neutrino and the recoil, respectively. In the frame where p'jn = 0, 
all v'^ - are proportional to Am^^^, so that no matter how exactly the average v'^ is defined, p'^jX" 
is second order in Amj^ and therefore negligible. Consequently, the phase reduces to the standard 
expression 

-E'l/'+p'l^^x" (15) 

that does not depend on the properties of the recoil. Due to Lorentz invariance, the phase must 
therefore be the standard one in any frame. 
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V. PLANE WAVES VS. WAVE PACKETS 

It is clear that, due to the Heisenberg uncertainty principle, neutrinos produced and detected 
in localized regions of spacetime must have a nonzero spread in energy and momentum. Many 
authors have argued that, to take this spread into account, a fully consistent theoretical treatment 



of neutrino oscillations requires wave packets 
wave packets are unnecessary 281]. 



Other authors have argued that 



Indeed, in most neutrino oscillation experiments, wave packet effects, such as wave packet 
separation due to different group velocities [9], are negligible, and simplified plane wave approaches 
can correctly predict experimental results. The deeper reason for this can be understood if we 
note that the coordinate space representation of a typical Gaussian neutrino wave packet, 

(X-V,t)2- 



V'j(x,t) oc e-^-^^o^+^P^o^'exp 



(16) 



is simply a plane wave, multiplied with an enveloping Gaussian. Here, the index j distinguishes 
different neutrino mass eigenstates, Pjo and Ejq are the average momentum around which the 
wave packet's momentum distribution is peaked and the associated energy; vj = pjo/Ejo is the 
group velocity, and ax is the width of the wave packet in coordinate space. The neutrino oscillation 
probability depends on the phase differences between wave packets associated with different mass 
eigenstates j. We observe that, at each fixed spacetime point, these phase differences depend only 
on Ejo and Pjo-^ They do not depend on ax, so they are independent of the precise shape of 
the wave packets, and, in particular, remain unchanged in the limit ax oo corresponding to 
plane waves. They do, however, vary over space and time. In the full wave packet picture, the 
enveloping Gaussian ensures that only space-time points along the trajectory x ~ vt, and located 
at the detector site, x ~ L, contribute to the oscillation probability.'^ In a plane wave approach, 
this has to be ensured by imposing x = vt and x = L by hand. 

Another argument not to invoke wave packets is the observation by Kiers, Nussinov, and 
Weiss [8] that a continuous flux of neutrino wave packets with identical momentum distributions 
cannot be distinguished from an ensemble of plane wave neutrinos, whose individual momenta fol- 
low the same momentum distribution. In fact, the density matrices describing the two ensembles 
are identical. 
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Note that interference can only occur between wave packet components located at the same point in spacetime. 
^ Here, v should be understood as an average of the individual Vj, which is a valid concept as long as the wave 
packets have not yet separated. Also note that |v| ~ c, with small corrections to this relation being a negligible 
second-order effect in the small neutrino masses. 
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FIG. 1: The Feynman diagram of neutrino oscillations. P and R are the parent particle and the recoil, 
respectively, while T is the target atom at the detector and S denotes the outgoing "signal" particles. 

Even though the above arguments show that plane wave approaches to neutrino oscillations can 
be justified, there are inconsistencies with the quantum field theoretic (QFT) formalism, where 
the production and detection processes are explicitly included in the calculation and the neutrino 
is treated as an internal line in a Feynman diagram (see Fig. [1]). Since energy and momentum 
are exactly conserved at the production and detection vertices, the intermediate neutrino can only 
be an on-shell energy-momentum eigenstate (as it has to be if it propagates over macroscopic 
distances) if its interaction partners — the external lines of the Feynman diagram — are also energy- 
momentum eigenstates. However, it is clear that the external energies and momenta in this case 
cannot be the same for diagrams involving different neutrino mass eigenstates. Since in quantum 
mechanics interference is only possible between amplitudes for processes that describe different 
paths leading to the same final state, this means that no oscillations are possible 
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A similar argument can also be invoked to show that, in the QFT treatment, entanglement 
between the neutrino and its interaction partners cannot persist asymptotically for t —?■ oo. Such 
asymptotic entanglement would lead to different external states in diagrams involving different 
neutrino mass eigenstates, so interference between these diagrams would be impossible. If we take 
the point of view that initially the neutrino is entangled with at least some of its interaction part- 
ners, we are led to the conclusion that these particles must be disentangled by further interactions 
for oscillations to take place. (This has been noticed previously in ref. [3].) Disentanglement is 
achieved when the particles interact with their environment. These interactions localize them and 
thus introduce momentum uncertainties that are usually large enough to allow external particle 
states entangled with different neutrino mass eigenstates to interfere in spite of their different 
energies and momenta. 
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With this in mind, one can justify a plane wave approach even in QFT. The procedure is 
to compute the amphtudes corresponding to different neutrino mass eigenstates and then sum 
them coherently (before squaring), keeping in mind that such coherent summation is only justified 
if the final state particles undergo further interactions that introduce energy and momentum 
uncertainties larger than the energy and momentum differences between the different entangled 
states. 

If one wishes to avoid such reasoning, one can also directly incorporate the effect of the localizing 
secondary interactions into the Feynman diagram computation by treating the external particles as 
wave packets with appropriate energy and momentum spreads. Different neutrino mass eigenstates 
will then simply couple to different (though usually overlapping) portions of the external wave 
packets' momentum distributions. Feynman diagrams involving different neutrino mass eigenstates 
can thus have identical external states, so that interference and therefore neutrino oscillations 
emerge. 

VI. CONCLUSIONS 

For the reasons we have explained, one must be cautious when following a plane wave approach 
to neutrino oscillation, whether or not the approach incorporates entanglement. A plane wave 
approach does not include the neutrino source and detector localizations that are physically essen- 
tial if oscillation with distance is to be observed. Consequently, such an approach cannot be used 
to discover all the implications of these localizations. However, the explicit inclusion of localiza- 
tion effects, which is automatic in wave packet treatments, can be replaced for some purposes by 
physical reasoning in a plane wave approach. Indeed, the oscillation wavelength derived from a 
plane wave analysis, with or without entanglement, can be perfectly correct. In particular, as we 
have shown, the wavelength that emerges from the analysis presented in this paper, which does 
invoke entanglement, is physically equivalent to the (correct) wavelength obtained by standard 
wave packet or plane wave treatments. Thus, one may continue to have confidence in the neu- 
trino mass-squared splittings Am^ that have been deduced by applying the standard wavelength 
expression to data. 
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